Black Holes 
The Theory 


Dr Dan Wilkins 
saint Mary’s University, Halifax, NS 


) 
SAINT MARYS 
v, UNIVERSITY SINCE 1802 
CRAQ Compact Objects Summer School, Montréal, August 2016 One University. One World. Yours. 


Overview 


Part 1: Developing General Relativity 
Part 2: Black holes and their properties 


Part 3: Accretion, from simple theory to messy 
reality 


Part 4: Using the theory to get close to black holes 
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Developing General 
Relativity 
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General Relativity - Why Bother’? 


e Newtonian gravity 
f= mg = —тУФ 
e Poissons equation 


V2 = AT Gp 


e No explicit time dependence — potential responds 
instantaneously to change in matter density over all space 


e Information travels faster than the speed of light! 
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Tne Equivalence Principle 


° Gravitational mass (coupling to gravitational field) 
vs. inertial mass (resistance to applied force) 


3 E d^3 ma > 
ana EO Ф — = —— VÈ 
J nav dt^ mq 


° Well-established experimental fact that ratio mg/m 
is constant for all particles 


° Trajectory of a particle through a gravitational field 
is the same, regaraless of the nature of the particle 
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Einsteins Elevator 


In a freely-falling (non-rotating) laboratory occupying a 
small region of spacetime, the laws of physics are those 
of special relativity — the (strong) equivalence hum 
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Special Relativity 


No preferred frame of reference 


Time passes differently for different observers 


AU = yAt Е 


Length contraction 


Ax 
Adz’ = — 
^y 
Mass/energy equivalence | = тос? 


Information cannot propagate taster than с 
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Spacetime Diagrams 


CI 
Future • теке 
а 
СА — Ат“ > 0 
МО 
Causal ü 
contact ° spacelike 
: E c At’ — Ax? < 0 
ж е Null (photon) 
Past c^ At? — Az? = 0 
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Lorentz Transformations 
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Lorentz Transformations 
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Lorentz Transformations 
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Gravity as Spacetime Curvature 


“Matter tells soace how to curve. Space tells matter 
how to move" — John Archibald Wheeler 


SAINT MARYS 
v, UNIVERSITY SINCE 1802 


One University. One World. Yours. 
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Gravity as Spacetime Curvature 
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Gravity as Spacetime Curvature 


“Matter tells soace how to curve. Space tells matter 
how to move’ — John Archibald Wheeler 
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The Equivalence Principle 
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The Equivalence Principle 
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Mathematical Toolbox — Vectors 


“ Vectors are geometric entities — exist Independent 
of any co-ordinate 

т? Сап be constructed trom 
basis vectors (components) 


у = 4)9 e, = „е, 
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Mathematical Toolbox — Vectors 


Vectors are geometric entities — exist independent 
of any co-ordinate 

Q2 а? Can be constructed from 
basis vectors (components) 


/ / 
v = ие, = u “e, 


A Obey transformation law 


la 
га __ Ox b 
U I ps U 
EA | Ол | 
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Mathematical Toolbox — Vectors 


Vectors are geometric entities — exist independent 
of any co-ordinate 
ШЕ u^ Can be constructed from 


basis vectors (components) 


/ / 
у = Vea = 0 е, 


A Obey transformation law 
L 


v^ = ail U 
т! i Әх? 


€1 " 
· Tensors are also geometric entities — generalisation ОГ, 
vectors to higher rank and transform as = = ЕУ = је 
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Mathematical Toolbox — Vectors 


Vectors are geometric entities — exist independent 
of any co-ordinate 
"EP Can be constructed from 


Summation convention: Where a dummy index 
appears twice (and only twice), once upper and once 
lower, it is summed over 
ues > uU va = У 1 
| €1 E КИ 
· Tensors аге also geometric entities – generalisation of 
vectors to higher rank and transform as = = ЕУ t E 
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Some Useful 4-Vectors 


Position x = (ct, x) = (ct, 2, y, 2) 
d Sa 
Velocity у = de (ci, JE (ye, Yu) 
dT 
Momentum р = Mov = (утс, ymv) = (Е/с,р) 
Wave vector p = hk k = (е Ju k) 


Current density j = ou j= (pee. j) 


Electromagnetic potential A = (ele, À) 


Proper time, т, is a scalar BE 
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Vectors in Curved Spacetimes 


Co-ordinates and Basis Vectors 
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Differentiating Vectors 


е LN gag 
| dr" Oxo ^ Ox? 


NV gue 


The covariant derivative: 


Ove 
Vav? — — + га yE 
Th on. a ге ag tl 
e connection: 5 = axe 
XL 
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Parallel Transport 


до“ | 
— =0inllat space 
du 


Vector parallel transported along a curve parameterised as x^(u) 
ду _ йө" " a ёа _ 4" -— „ Ова 41° 
du ди“ Ë du da Әх du 


dv _ (dut | apa ARE. Dv? 

в d uu вс? т | Ca Tu O for parallel transport 
№ 

Dv° EE 
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Generalise the flat space line element ds? = (dx)? + (dz?)? 
ds? = guu (82 )- + goo (dz)? + 2901 (dz!) (da^) 
In general ds? = gapdx*dx” 


The metric tensor, gap, encodes all of the 
geometry of the spacetime, with 6а: ep = дар 
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Contravariant vs. Covariant 


• For any set of basis vectors {ea}, can define the 
dual basis (e^) such that e, : eb = ô? 


a 
| | | У = 79 e, 
° Can write a vector м іп terms of either. oa 


Va e 
° o^ are the contravariant components of v 
° ua are the covariant components of v 


° Since eu : = дар, Va = gay v" and 09 = gw, 


9°") = [ge] ` pes 
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Length of a Vector 


е. The line element is the length of the (infinitesimal) 
displacement vector 


аз? = ga dax? dx” 
° Apply to other vectors 
КЕ = V- V = gay v = VV, 


• Dot products are invariant (same in all frames/co- 
ordinate systems) 
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Making an Observation 


e.g. 4-momentum p which has 
components in Minkowski 
Space: 


ір” I (Ep) 
Observer measures energy 


Е = реу=р и 


• Given a vector, у, an observer measures the components against their 
own set of basis vectors (flat, Minkowski space in their rest frame) 


° They are at rest in their own frame, so their timelike basis vector must be 
parallel to their 4-velocity 


• Can construct an orthogonal set of vectors around this to give x, m А 
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Geodesics 


= B 
e Shortest path between 2 points s = J ga i? du 
A 


* Tangent vector (the momentum) is parallel- 
transported along the geodesic 
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What Is the appropriate 
Lorentz transformation for the 
density of a cube of matter? 
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Describing me v апе 
KO 


ПЕ 


° Density increases by factor y —transforms as a rank 2 tensor 
TU = poutu” 
° The energy-momentum tensor for a perfect fluid 
Т!” = (po + р/с”) utu” — pg” 
e To conserve energy and momentum У„ТР” = 0 
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Describing the Curvature 


e |n flat space, it doesn't matter which order you perform a 
covariant derivative (because in flat space these are just 
partial derivatives) УМ „Уро“ = VV и“ 


• Can measure the curvature of space by changing the 
order of covariant derivatives 


Ve Va — МУ Va = В‘ ла 


• Ra 15 the rank 4 Riemann tensor (or curvature tensor) 


• Can also contract it (sum over indices) to give the Ricci 
tensor (rank 2) and scalar 


__ рс _ pa 
Rab = abe R = На ASS 
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Tne Einstein Field Equations 


° Curvature of soacetime is related to the matter 
content and energy-momentum tensor Is rank 2 


Gag RE 


° Needs to reduce to the Newtonian limit (terms no 
higher than linear in second derivatives of metric) 


° | is symmetric so G must be symmetric 


1 SiG 
Ка — SOA = A 
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Tne Einstein Field Equations 


° Curvature of soacetime is related to the matter 
content and energy-momentum tensor Is rank 2 


Gag RE 


° Needs to reduce to the Newtonian limit (terms no 
higher than linear in second derivatives of metric) 


° | is symmetric so G must be symmetric 


1 SiG 
Ка — SOA + Agw = —— in 


SAINT MARY'S 
GO UNIVERSITY SNC 1802 


Е! П Stel na 1 9 1 6 One University. One World. Yours. 


One University. One World. Yours. 


Gravity around a Point Mass 


The Schwarzschild Metric | 


9 —1 
ds* = c° (1 — +) dt? — (1 — +) dr? — т? (40° + sin? дар) 
т T 
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Gravitational Redshift 


In an observer's rest frame, Spacetime is described 
by the Minkowski metric (using their co-ordinates) 


Their own worldline ds? = dr? 
2 
сіт? = c° f = +) dt? 
T 


For observers closer to the black hole, proper time 
elapses more slowly 


Photons climbing out of the potential are redshifted 
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Geodesic Equations 


In the equatorial plane of the Schwarzschild geometry (Ө = x/2) 


1-28) г (1) 


с2 (1 = >) i? — (1 u У B $20 2? _ c*, massive particles 
d 4 0, photons (2) 


"“ф-һ (3) 
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Constants of Motion 


r 


h is identified with the angular momentum 
• To understand k, consider the 4-momentum of the particle 


e An observer with 4-velocity u measures the components of the 4-momentum 


as [p^] = (E, p) 
e Energy is measured as, № = р · u which for an observer at rest at infinity 


Е = дор” = ро = ke? 


e kis Ше identified with the energy of the particle in its orbit k = 
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Energy and Shape Equations 


* Energy equation — substitute (1) and (3) into (2) 


p t = GUK — 1) 


h? f ú = |. 2GM 


сет г 


• Shape equation — substitute (3) into above and 
differentiate (u = 1 / r) 


Фи GM зам, 
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Radial Motion 


Setting 40 = аф = 0 
energy equation becomes 
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Radial Motion 


Setting 40 = аф = 0 
energy equation becomes 
GM 


r2 


į? = @ (k? — 1) + 
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Radial Motion 


Setting 40 = ар = 0 LA 
energy equation becomes “ d, pups ЕЗИ 3 


1—2 
E 


t drd 


r 


і — со as r— 2и 


Ку hn 


| | | | | | | | | | | | 
4 6 & 10 
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Т you were to drop objects 
Into a black hole, what would 
it look like some time later’? 
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Orbits 


° Set dr = d0 = 0 for orbits in the equatorial plane 


O 2GM\ 2GM a-a 
e Energy equation = ( - “2, )- ~ = с(К — 1) 
| М М 
° Shape equation u= m = и? 
np ET 
г — Зу (1—3u/r)? 


do ЕМ 1 
dt r3 ES SAINT MARY'S 
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Stability of Orbits 


Newtonian orbits 
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stability of Orbits 


2 с^ ; 
zer Ede 5 (0-1) 


In Schwarzschild 
Spacetime 


| Innermost stable circular 
orbit ISCO) at r = би 
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(In simple terms) what would 
an accretion disc look like 
around a black hole”? 
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Motion of Photons 


с^ (1 — = i? — (1 __ > B pe — r2 22 B C^, massive particles 
Џ 4 0, photons 


° Energy equation becomes 


° And the shape equation 
du | зам, 
аф? с? 
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Gravitational Lensing 


From the shape equation for a small perturbation Au 
to the straight line path 


Pu E = 
|22 Tu с^ i | | | 
ET ER SAINT MARYS 
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Gravitational -ensing 


D 
Ве = ADrs b= Вр 

Ds 
pa СМ p, Pr ө. НЕ _ 
i с^ > E Dg 


IGM Dis 
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Gravitational -ensing, 


] / 
nS ARE r) 


ћаћ 


о? +2 (T 
~ dp N 


о +14 | Re ey sane MANG 
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Magnification = 


Light таме! Delays 


R Рә 


(R — 2u)^ R^ 


~200us for echoes between Earth and Venus T SAN MANG 
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Circular Motion of Photons 


d'u 3GM о 
dio? и = и 


| | 3GM 
2 circular orbit at r = 5 

| С 
Determine stability trom Ver | ||| || ||| 


energy equation 


Уен 


Maximum in potential 
e Unstable orbit 


6 7 8 
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Tne Equivalence Principle 
(revisited) 


2 2 Detector 
At what frequency 


E is the photon 
detected? 
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Tne Event Horizon 


= 


| 2 
ds? — с^ __ SE dt? = = aie dr? — „2 (40? 4 РЕ 949) 
| T 
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Rotating Black Holes 


The Kerr Metric 


2 4 0 p? 2ua?r sin? 
ds? = (2 (1 - ZZ dt? + — 0 ар = dr? — paf? — G ЕЕЕ 3 sin? даф? 
p 


GM 
Y 


и = 
p^ = т? + а? cos? 0 


A = г? — рт + а? 
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Rotating Black Holes 


The Kerr Metric 


ds? = Ky Фар -— — de GO 042)? Е аг ú Рав 
GM 
и= => 


A = r^ —2pr + a° 


= (” + а?) — a° Asin? 0 
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Rotating Black Holes 


The Kerr Metric 


2A 572 я 20 А 
GM 
— E 


р = r? + а? cos? 0 
A — r^ -2иг фа 
Y = (r? а?) — а2А sin? 0 | 


2ar 


Ш = Уу? 
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Rotating Black Holes 


The Kerr Metric 


2 A Y? sin? 0 2 
GM 
р =r? + а? cos? 0 
е J = Мас 
5? = (nee) — PA sin? af 
2ar 
uno] 
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Structure of a Kerr Black Hole 
Event horizons ~ 
rp = и + (u^ = а?) | 
_ = и (м? — a?) 


м мін 


+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 
Ф 
Ф 
+ 
+ 
+ 
Ф 
+ 
+ 
+ 
+ 
+ 
4? 
+ 


ча 
а 
ча 
а 
ча 
"aa, 
.. 
ча 
"a 
.. 
.. 


Ergosphere 


Ee 
"a 
ча 
.. 
а 
а 
tea, 
а 
ча 
à 
.. 
ча 
а 


Ring singularity о \ 
Ga | Pot = И ER (ГЕ GO д2 cos’ 0) 


rs- = и — (ГЕ — a? cos? 0) 
Stationary limit surfaces 
SAINT MARY 5 


— (Infinite redshift) Se isse = 


stability of Orbits 


8 
Retrograde, a < 0 
6 
4 Prograde, a > 0 
2 


02 04 0.6 0.8 10 ° 


One University. One World. Yours. 


What would you see as 
you fall into a black hole? 
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Accretion 


Powering some of the brightest objects in the 


EIZ 


Universe ANTA 
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Accretion Discs 


To transfer matter inwards, need to transfer angular 
momentum outwards 


GM 


„тетш e 7- 120 =VGMr 
Г 
Total angular momentum transfer rate J = MVGMr 


If some fraction 6 can be accreted from inner edge 


j= м (VGMr -вусмп,) 


Shakura & Sunyaev 1974, Pringle 1981 


Accretion Discs 


12h 


Angular momentum in transferred by viscous stresses 


мета 
Rate of change of angular momentum equal to torque 


s,(2rr2h)r = М [v GMr — Ву GMrin) 


M(VGMr — BVGMrin) 


Sp — Алт? В, 


Shakura & Sunyaev 1974, Pringle 1981 
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Disc Emission 


Work Done = Force x Velocity 


Power _ Force x Velocity ` Force i Velocity 
Unit Volume . Volume |J Area Radius 
| du; d) 3 du (т) _ dQ 
Q= $e dr pA dr 9 "e dr dr ци dr 
So the power output per unit volume of the disc 
к d M(VGMr — BVGM Trin) 
С = с EEE 
nr<h 


Shakura & Sunyaev 1974, Pringle 1981 
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Disc EmIssion 


° From the power output per unit volume, can 
integrate through the height of the disc to get the 
flux (power per unit area) 


° Half comes out the top of the disc and half out the 


bottom 

1 . 

F = Әһ 
> Q 

3GMM m 
> 1 — 
dud (1 p GA 
ESU As 


Shakura & Sunyaev 1974, Pringle 1981 One University. One World. Yours. 


The Disc spectrum 


Assume disc is optically thick so each part of the disc 
radiates as a black body 


FE 


f 


7/1 | __ГЕ(т) 1 
LY, те) = |=) 
PA 
b 


/ 
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Tne Viscosity 


° Viscosity expressed in terms of the hydrostatic 
pressure in the disc sí = op With а < I 


° a 0.1 to produce accretion rates high enough to 
produce observed luminosity 


• Cannot be explained by viscous friction... 
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Accretion Efficiency 


Assume the gravitational potential lost between 
Infinity and the innermost stable orbit is radiated 
away (all radiation from within is lost into black hole) 


Newtonian orbits: 


p СММ 1 > GMM 
T Z 2r 
€ = : r= би 2 
М с? 12 
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Accretion Efficiency 


Assume the gravitational potential lost between 
Infinity and the innermost stable orbit is radiated 
away (all radiation from within is lost into black hole) 


Schwarzschild spacetime: 


E — >Ñ E 
Е = ССО k = e=] - [ 
тос? тос? 
k — d= Зрјт) € = 5.7% for a Schwarzschild 
(1 — Зи/т)г black hole 
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How does the efficiency 
of a rotating black hole 
compare"? 
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Radiatively Inefficient Accretion 


| log R (Rs) ° 
Narayan 1998, Narayan & Quataert 2005 
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Radiatively Inefficient Accretion 


! log R (Rs) ° 
Narayan 1998, Narayan & Quataert 2005 


Building a more realistic picture of 
the accretion disc 


^ eau. "ow 
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Magnetohydrodynamics 


If the accretion flow is made of an ionised plasma, it can be 
described by magnetohydrodynamics (MHD) 


e Continuity equation (mass conservation) 


Op 
eda JA. — 0 
• Charge conservation 
Og — — 
—+V-j=0 
TREE 


• Momentum equation 


Où 4 4.2 - 
(Жа voit, 
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Magnetohydrodynamics 


Electric and magnetic fields obey Maxwell's equations 


э.б 4 
60 
a. 2 OB 
VAE = 
Ot 
У.В=0 
> a ‚ LOE 
A В = — 
M TIT O 


And Ohm's law j = e(E + л B) 
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Ideal Magnetohydrodynamics 


° Charge neutrality — charge imbalance causes motion 
of electrons to balance it out. Charge imbalance 
oscillates at the plasma frequency - effective 
shielding length 


mE 1 8E "m т 
е Non-relativistic flow 2 OF <VAB V ^ В = uoj 


° Infinite conductivity — can neglect diffusion of the 


magnetic field ОВ Е = 1 - 
lid = У A (t A В) + —V*B 
дъ под 
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Ideal Magnetohydrodynamic 
Equations 


° Taking the fluid equations and Maxwell's equations with the 
assumptions of ideal magnetohydrodynamics and eliminate E 


СР +. (pil) = 
ДЄ i Vi) =JAB- Sp 
Falan B) = 22 
V ^ B = и] 
p= Kp! у= 5 е 
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The Magnetic Field 


V ^ B= no) 
fing ec = —]+(B-V)B 
а= |9 (5) 08-9 
В? 
e Magnetic pressure Pmag = -- 
2 uo 


e Magnetic tension from (B-V)B 
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MHD Waves 
— EEE | 


Longitudinal 
Up = Cs ШШ NM DU DU › 
— 
Along field lines 
Transverse 
Up = VA 
— 
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Magnetic Flux Freezing 


Flux through a surface co-moving with the fluid (velocity ш) 

4Ф е 4 a 4 

tes | Бе an -48 | B(t) -dS 

dt dino m | S(t) с) | 
The flux through the total closed surface (£ is 
the edge, dx a line within it) must be zero 


|| B(t + dt) - 8+ | В (t + dt) - gura: | B(t + dt) d$ = 0 
S (t--dt) S(t) 


Use Stokes' theorem to integrate around the curve dx that encloses 
the surface S(t) and use MHD equations 


dd > do 1 > – 
E - | ‚459 = = | сој . 45 
dt S(t) dt S(t) N Ho? 
Em SAINT MARY'S 


AAA dé | | я 
If conductivity is infinite — = 0, flux moves with fluid ҚУ) asin 
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Producing an X-ray Continuum 


Comptonisation | 
Continuum 


Thermal distribution of particles 
Temperature T = Multiple (inverse) Compton 


on from disc 


Average energy change of photon per Compton scattering 


(=) _ АКТе — hu 


€ тес? 


The Compton optical depth y = 


T. aa 
2 SII SAINT MARYS 
-z шах (27 г. Қ SAINT M \КҮ 5 


e One University. One World. Yours. 
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State of the Art: GRMHD 


° Need a covariant description of the electromagnetic field 
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GRMHD 
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simulation Diagnostics 
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Summary 


• Space time around black holes described by general relativity 
• Particle orbits, passage of light, redshifts 


• Accretion physics describes the inflow of gas and how luminous 
objects are powered 


° Simple models, MHD, GRMHD 
* Theoretical models let us interpret observations 
* How are the most luminous objects in the Universe powered? 


How does the corona work? How does it vary? How is it 
governed over cosmic time? 
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